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Abstract 



In this paper we study long distance modifications of gravity obtained by 
considering actions that are singular in the limit of vanishing curvature. In par- 
ticular, we showed in a previous publication that models that include inverse 
powers of curvature invariants that diverge for r ^ in the Schwarzschild ge- 
ometry, recover an acceptable weak field limit at short distances from sources. 
We study then the linearisation of generic actions of the form L = F[R, P, Q] 
where P = Rf.uR'"' and Q = R^^XaR^"^" ■ We show that for the case in which 
F[R, P, Q] = F[R, Q — 4P], the theory is ghost free. Assuming this is the case, in 
the models that can explain the acceleration of the Universe without recourse to 
Dark Energy there is still an extra scalar field in the spectrum besides the mass- 
less spin two graviton. The mass of this extra excitation is of the order of the 
Hubble scale in vacuum. We nevertheless recover Einstein gravity at short dis- 
tances because the mass of this scalar field depends on the background in such a 
way that it effectively decouples when one gets close to any source. Remarkably, 
for the values of the parameters necessary to explain the cosmic acceleration the 
induced modifications of gravity are suppressed at the Solar System level but can 
be important for systems like a galaxy. 
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1 Introduction 



Theories that can accommodate a long distance modification of gravity are interesting 
phenomenologically, but also from a purely theoretical perspective. Phenomenologi- 
cally, because although Einstein gravity has passed many tests in the Solar System, 
the issues of Dark Matter and Dark Energy might well have something to do with a 
non-standard behaviour of gravity at large distances. And from a theoretical stand- 
point any such theory is worth consideration because it has proved very hard to find 
consistent, generally covariant theories of this type. In this article we report in what 
we believe represents a consistent, large distance modification of gravity that, besides 
being able to explain the acceleration of the Universe, can also have interesting impli- 
cations for the Dark Matter problem. These modifications of gravity are generated by 
including inverse powers of the curvature in the action, or, more generally, considering 
actions that are singular in the limit of vanishing curvature. These models were first 
proposed to explain the cosmic acceleration without the need for Dark Energy [Ij, but 
it was soon realised that the simplest possibility {AC oc 1/R) can give problems with 
the weak field limit j2j. This can be easily understood noticing that models whose 
Lagrangian is an arbitrary function of the scalar curvature are generically equivalent 
to scalar-tensor theories of gravity. In the models that can explain the acceleration of 
the Universe by introducing inverse powers of the scalar curvature the corresponding 
scalar has a mass that is proportional to some positive power of the scalar curvature. 
One can expect then that in regions where the scalar curvature is very low this excita- 
tion generates a new force that would be in confiict with experiments. And although a 
rigorous quantitative computation of the predictions of these theories for Solar System 
measurements is still lacking in the literature it is not clear that they constitute a 
viable alternative to Dark Energy^. 

But the situation is quite different if we consider inverse powers of other invariants. 
We showed in a recent publication [3] that one can generically recover an acceptable 
weak field limit in these models at short distances provided one includes inverse powers 
of curvature invariants that do not vanish for the Schwarzschild geometry. In particular, 

-•^For instance, one of the conseqences of this scalar field is that one would infer a different value of 
Newton's constant from Newton's potential and light bending. 
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we considered the action 

S = I d'^Xy^^ 



with positive n and where 



,,4n+2 

R- ^ 



[aR^ + bP + cQY 



P = R^,R^^ and Q = R^^xpR"""^' ■ (2) 

Here Gat is Newton's constant and /i a parameter with dimensions of mass. In this 
theory corrections to standard cosmology will only occur at the present epoch if we 
take /X ~ i^o and it has been shown that it has cosmologically interesting accelerating 
solutions [4^ that provide a good fit to the SN data |5_^ without recurring to Dark 
Energy. But we also showed in [3] that beyond a source dependent distance given by 



(3) 



where M is the mass of the source, gravity is modified. Moreover, we computed the 
first correction to Newton's potential in an expansion in powers of r/r^ yielding 

12n+8^ 



V{r) 




^, (4) 



where a ~ 10^^. It is remarkable that with the values of the parameters necessary to 
explain the cosmic acceleration, this characteristic distance puts the departure from 
Newtonian gravity well under control at Solar System scales, but can make it important 
for larger systems like a galaxy, hence the potential implications for the Dark Matter 
problem. But the theory under consideration involves up to fourth order derivatives of 
the metric. The existence of these higher derivatives suggests that one would always 
find ghosts in a linearised analysis. In fact, it can be argued by considering the Cauchy 
problem and the gauge symmetries and constraints of the theory (see for instance [6 ) 
that general actions involving arbitrary powers of the Riemann tensor will contain at 
most eight degrees of freedom^: two in the usual massless graviton, one in a scalar field 
and five in a ghost-like massive spin two excitation. The problem with the ghosts is 
that one has to accept negative energy states in the theory or lose unitarity [7j. In the 
former case the background in which they appear is expected to be unstable at least 



^If we also included derivatives of the Riemann tensor in the action even more degrees of freedom 
would show up [H]. 
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towards the generation of small scale inhomogeneities (positive and negative energy 
particles) since that configuration would dominate the phase space volume available 
to the system. Therefore the consistency of these theories has been put into question 
|8|9j . In this paper we will show that with a suitable choice of parameters {b = —4c), 
the theory obtained from the action is ghost-free. In fact we will show that generic 
Lagrangians of the form C = F[R, Q — 4P] are ghost free. As expected, in this case 
we just add a single scalar degree of freedom to the gravitational sector, on top of the 
massless spin two graviton. For theories of the type that can explain the acceleration 
of the Universe without Dark Energy, the mass of this excitation is ~ H. Nevertheless 
we still recover an acceptable Newtonian limit at short distances from sources because 
the linearisation of the theory breaks down in this domain, and the extra scalar degree 
of freedom decouples. 

At this point it is convenient to recall the properties of other theories that offer 
long-distance modifications of gravity, since we will see that all of them share common 
features that are expected on general grounds for this kind of theories (see for instance 
the discussion in ^^). The simplest way to modify gravity at large distances would 
seem to be to give the graviton a tiny mass, rrig. This is not as simple as it seems, 
because when doing so we also introduce more degrees of freedom in the gravitational 
sector, and this has some unexpected consequences. If we study the theory at a lin- 
earised level, the only ghost-free mass term we can add is the Fierz-Pauli mass term 
When adding this mass term one finds out that the linearised theory does not 
recover the results of Einstein gravity when we take the massless limit 0. This 

is the so-called vDVZ discontinuity and although it was first interpreted as a 
failure of massive theories of gravity, it was later interpreted as a failure of the lineari- 
sation |13|14j . This is so because the new extra degrees of freedom that the Fierz-Pauli 
term introduces have couplings that are singular in the massless limit. This can be 
qualitatively understood noticing that the "kinetic term" of the extra degrees of free- 
dom comes, roughly, from the Fierz-Pauli term. Once the fields are rescaled according 
to their canonical normalisation, the mass of the graviton appears suppressing non- 
renormalisable operators in the expanded action. For a small graviton mass one finds 
that the "cut-off" of the theory is very low and the linearisation breaks down when 
one gets close to a mass. Therefore the results obtained from the linearised Lagrangian 
are not valid at short distances in the spacetime of a spherically symmetric mass. The 
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distance at which the hnearisation breaks down is usually called the Vainshtein radius 
since he was the first to argue that in a fully non-linear generally covariant theory that 
reduces to massive gravity upon linearisation one should be able to check that the vDVZ 
discontinuity is absent J^l- Continuity in the massless graviton limit should then be 
achieved non-perturbatively. Unfortunately it is not possible to obtain this limit if we 
don't have a full non-linear generally covariant completion of the Pauli-Fierz term, and 
the issue of the consistency of massive gravity remains open. 

The story is very similar in other more elaborate, generally covariant models that 
also yield a long-distance modification of gravity. For instance models involving Lorentz 
violating condensates for some extra fields have been considered [T5] to achieve this 
kind of modifications, and these models have been shown to be consistent. Also theories 
obtained from braneworlds like the GRS JH] and DGP J7j models show this kind of 
phenomenology. But while the GRS model presents ghosts upon linearisation, these are 
not present in the DGP case. This last model is probably the best studied long distance 
modification of gravity so far. It is related to massive gravity, since the graviton 
acquires a "soft" mass in its propagator. In this model the strong coupling scale and 
the nonlinear effects that imply the absence of the apparent vDVZ discontinuity have 
been studied in detail (see e.g. jlH] and references therein). The situation can be 
described as a shielding of the would-be strongly coupled modes by non-linear effects 
at short distances. In this case gravity becomes higher dimensional at long distances. 

As we have mentioned, the model we will deal with shares many common features 
with these theories: looking at the linearised action on de Sitter background there is an 
apparent discontinuity in the would-be Einstein gravity limit (/i 0) due to an extra 
scalar degree of freedom that becomes massless and would seem to make the model 
not realistic. That is again resolved once one takes into account that this linearised 
version of the theory is no longer valid in regions where the curvature is large. The 
extra scalar degree of freedom gets shielded and decouples at short distances, affecting 
the dynamics only at long distances where the background curvature is small. Our 
model, however, does not produce a massive graviton: the spin two excitations of 
the metric remain massless. In the next section we will briefly summarise how and 
why do we recover the Schwarzschild solution at short distances, r -C r^, making 
explicit the shielding mechanism. In the third section we will discuss the linearised 
action for generic Lagrangians of the form £ = P, Q] and we will show that when 
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F[R,P,Q] = F[R,Q — 4P], the theory is ghost free. In the fourth section we will 
clarify the range of validity of the linearised action, having in mind theories of the 
type (0) that can explain the acceleration of the universe without Dark Energy. We 
will see that the "cut-off", or strong coupling scale of the linearisation in vacuum is 
given by ~ {MpH^Y^^ and we will obtain the Vainshtein radius for an object of 
mass M as ry ~ (G^M/ H^Y^'^. This does not mean that the theory loses predictive 
power at higher energies or shorter distances, but one has to use non-perturbative 
methods or a different kind of expansion to obtain the results. Section five contains 
the conclusions. In the Appendix we offer the technical details of our derivations: we 
give a gauge invariant decomposition of the degrees of freedom contained in the metric 
as well as the action expanded up to bilinear order in terms of these excitations and 
their propagators. 

2 Short distance solution 

In this section we will briefly summarise the results we obtained in pj. As we said we 
started with the action (^. To study the Newtonian limit we want to find spherically 
symmetric solutions to the equations 

G^, + /i^"+2^^, = 0, (5) 

where G^j^ is the usual Einstein tensor and ^^^^"^H^^ is the extra term generated by the 
inverse powers of the invariants in the action. The main point is that, when evaluated 
in the Schwarzschild solution, the extra term that appears in the equations of motion 
is of order 

^4n+2^.(o) _ [GM/r^){rlr,f^+\ (6) 

where was defined in the introduction. This shows that also G^j, tends to zero as 

r — > and we indeed approach the Schwarzschild geometry in this limit'^. At small 

distances, r ^ r^, we can then consider a small perturbation of the Schwarzschild 

solution and solve at first order in the perturbations, obtaining the correction to the 

Newtonian potential presented before. This amounts to an expansion of the exact 

solution in powers of r/r^, and is therefore only valid in the r <^ region. Notice that 

■^Remember that the Schwarzschild solution is the only spherically symmetric solution to the equa- 
tion Gfj,i, = 0. 
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an expansion in powers oirjrc signals a breakdown of the perturbative series in powers 
of Gjq. This feature is expected in any theory yielding an infrared modification of 
gravity ^lOj and it is related to the fact that the extra degrees of freedom introduced by 
the modification are shielded at short distances. The reason of the shielding behaviour 
is clear in this case: the Kretschmann scalar Q, when evaluated in the Schwarzschild 
solution reads Q = , while both P and R vanish. So, assuming c is not zero 

in (PJ), the term ^'^^^'^H^y has a shielding prefactor for spherically symmetric solutions 
that goes like ~ (r/rc)^""*"^, effectively switching off at small distances the modifications 
of the solution induced by this term. So we see that gravity is modified only at large 
distances in this theory for the same reasons why cosmology is modified only at late 
times, when the mean curvature of the universe falls down to a value i? ~ yU^. 

However, although we have seen that this modification of gravity carries with it a 
shielding mechanism for massive sources, we have not seen what is it shielding us from. 
One important requirement for any proposed modification of General Relativity is the 
existence of stable phenomenologically viable cosmological backgrounds that are free 
of ghosts and tachyons, or where, at least, their effects can be reduced to unobservable 
levels locally. For investigating this issue it is convenient to study the linearisation of 
the theory over one such background. This is the goal of the next section. 

3 Linearisation of modified gravity 

In this section we are going to study the linearisation of a generic modified theory of 
gravity. In the case of actions of the type ((T)), this expansion can be thought of as 
complementary to the expansion over a Schwarzschild background that we considered 
previously, one being valid at short distances and the other only at large distances. 
But as we said let us consider here a generic action built with the scalars i?, P and Q: 



The background over which we consider the linearisation should be a solution to the 
equations of motion in vacuum. If we consider at the moment constant curvature 
maximally symmetric spacetimes the possible values of the curvature can be found as 
the solutions of 




(7) 
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where Fq = OqF, etc... 

To determine the stabihty of these solutions and the particle content of the theory 
we should expand the action in these spacetimes up to second order in the fluctuations. 
It can be seen |6|9j that this expansion will be the same as that obtained from 



S = id'^x^/^- 



(9) 



where C^yXa is the Weyl tensor and we have defined 

A = ^F-RFr + i?2 (FnR/2 - Fp/A - Fq/6) + i?^ {Frp/2 + Frq/3) 

+R' {Fpp/8 + Fqq/18 + Fpq/6))^ (10) 
5 = {Fr-RFrr-R\Fpp + 2Frq/?>) 

-R' {Fpp/4 + Fqq/9 + Fpq/3)>„ (11) 
mo 2 = {{3Frr + 2Fp + 2Fq) + R{3Frp + 2Frq) 

+R' (3Fpp/4 + Fqq/3 + Fpq))^ (12) 
m^' ^ -{Fp + AFq)^. (13) 

Here < ... >o denotes the value of the corresponding quantity on the background. So 

we see that the situation for the perturbations here is the same as in conventional 

Einstein gravity supplemented with curvature squared terms. It is well known [6^ that 

in this case the gravitational sector contains a massive scalar excitation of mass ~ mg 

and a ghost-like massive spin two field with mass ~ m^, besides the usual massless 

graviton. So the ghost will be absent and there will be no ghost-like unstabilities in 

the theory as long as ^ = 0. In particular it is evident that actions of the form 

F[R,P,Q] = F[R,Q — 4P] are ghost free (see also [E]). But this expansion is only 

valid in constant curvature spacetimes. We would like to consider expansions in other 

FRW backgrounds that are also solutions of the equations of motion in vacuum, since 

for some phenomenologically interesting models de Sitter space is unstable (the scalar 

is a tachyon) and the late time solution corresponds to a power-law FRW cosmology 

jlj. We show in the Appendix that the absence of the ghost for actions of the form 

F[R, P, Q] = F[R, Q—iP] also translates into these backgrounds. This is not surprising 

since, locally (at scales much less than H~^), there should be no significant differences 

between these spacetimes^. Moreover, we give in the Appendix the decomposition of 

^One way to think about the absence of the ghost is that its mass has been taken to infinity. Its 
absence is then an "ultraviolet statement" that should not depend on the details of the background. 
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the metric and the hnearised action in terms of the physical gauge invariant excitations 
and their propagators. Using then the equation of motion for the extra scalar excitation 
it is straightforward to address the stability of the de Sitter background. This scalar will 
satisfy the Klein-Gordon equation with a mass given by (see eq. ()86|) of the Appendix) 



For perturbations with zero spatial momentum we then have: = 00^^*™''*- Taking 
into account the definitions, eqs. ()65l74|) . this results in a spacetime metric 



From this we can infer that the background de Sitter metric is stable for > —9H'^/4, 
consistent with the results of for F[R) theories. 

4 The limits of linearisation 

We will turn now our attention back to the original action (^. We have seen in the 
previous section that in the case b = —4c the theory is ghost free. But looking at the 
excitations over the vacuum it would seem that even if there is no ghost, the theory 
under consideration, containing an extra scalar degree of freedom with a tiny mass, 
would be ruled out by fifth-force experiments and Solar System observations, as was 
concluded in |9| . We will see that this is not in fact a problem, if one takes into account 
that the linearised action we are using to extract these conclusions has a very limited 
range of validity. In the expanded action higher order terms will be suppressed by 
inverse powers of the background curvature. This means that this expansion of the 
theory will break down when the curvature of the fiuctuations is extremely small, of the 
order of the background one. In regions of high curvature we can not trust the linearised 
action obtained from Q. So for a spherically symmetric mass there is a radius such 
that inside it we can no longer trust the results obtained from our expanded action. 
This is analogous to the Vainshtein radius of massive gravity. In the next subsection 
we will argue that the strong coupling scale in vacuum for theories of the type (Q) 
(with yU ~ H), or the cut-off of its linearised version, is A, ~ {MpH'^f^, where 
Mp = {IQtxGn)^^. This strong coupling scale corresponds to the scale suppressing 
non-renormalisable operators for the canonically normalised excitations. Notice that 




(14) 



8 



this does not mean that the theory loses predictive power above this scale. We know 
the ultra-violet completion of the theory, but at energies above this scale we should use 
non-perturbative methods or a different kind of expansion. For spherically symmetric 
solutions this gives a Vainshtein radius given by ry ~ {GNM/H^f^. In subsection 
4.2 we will explicitly check this result by computing the higher order corrections in 
spherically symmetric solutions, and we will find that higher order corrections do indeed 
become dominant at distances less than ry. 

4.1 Higher order corrections: the strong coupHng scale and 
the Vainshtein radius 

To obtain the strong coupling scale, we have to figure out the scales of the coefficients 
that go in front of the higher order non-renormalisable (marginal) operators, after 
canonical normalisation of the fields. We will assume that our Lagrangian is given by a 
function F[R, Q—AP\ and we are in de Sitter spacetime for definiteness, with R = 12H^. 
It will be convenient to write F[R, Q - 4P] = F[R, U], where U = 5i?V6 + Q - 
because U is zero on the background. We consider as usual the weak field expansion 
in powers of h^i, where g/^i, — 51°^, + h^i, and the background metric is 

ds"^ = glJx^dx" = -dt^ + a{tf{dx^ + dy"^ + dz^) (16) 

with a{t) = e^*. Retaining only the highest dimension operators at every order (these 
will dominate at high energies^), the expansion in /i^,^ of U roughly looks like: 

u ~ {d^hf + {{d^hfh + d^hdhdh) + {{d^hfh^ + ...) + .... (17) 

We have suppressed all the indices so we write d to denote any space or time derivative 
and h for any h^i, . For the expansion of the Ricci Scalar we have: 

R ~ 12i/2 + a(t)-='/' (^(V' + 3i/ao + yi/')</.+ (3ao2-2V' + 3//9o-^i/>^ 
+ {hd'^h + dhdh) + {{d^h)h^ + {dhfh) + . . . , (18) 

where is the three-dimensional Laplacian and = a(t)-2V^ We write the full first 
order term explicitly (using the gauge invariant decomposition given in the Appendix), 
because the scaling dimension of this term is not what it seems. At first sight this 
^Remember that by "high energies/momenta" we mean here H. 
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term seems to go like d^h at high energies. However, as we show in the Appendix the 
propagators are such that when doing perturbative calculations, it is more convenient 
to replace (j) with cf) (see eq. fl97p in the Appendix). This replacement diagonalises the 
propagators and cancels out the d^r and the Or terms in the first order term. From 
eqs. ()97ll00p we find now for high energies/ moment a: 

i?«^a(t)-3/2(^V20+(...)ifV) . (19) 

Here we have considered r and separately, since, as we will see later, the scaling of 
the (f) propagator is non-canonical. 

Using the general forms of the expansions of U and R we can now write the general 
form of the higher order terms in the expansion of the action, again only retaining the 
highest dimension operator at each order. This expansion will take the form 







+ {Fuur}^ {d^hYR^'^ + ( Fuuu )^ {d'hf + ...), (20) 

where S**^^-* is the action expanded up to second order in h (see the Appendix). In 
order to estimate the strong coupling scale, we still have to canonically normalise the 
fields. We show in the Appendix that for high energies uj and high momenta fc, the 
propagators of the six gauge invariant modes contained in h^y, go like: 

1 

< XijXij >o 

< WiWi >o 



M2(u;2 _ p) 
1 



< rr >o ~ (21) 

M2(cu2 - A;2) ' 



< 00 >o 



H 



2 



M2fc4 



One can read off now the canonical normalisation of the excitations: h = h'^/Mp, 
except for 



So in general, at high energies/momenta, the contribution of in a certain higher 
order term will be subleading with respect to the contribution of the other excitations. 
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However, as we have shown exphcitly, there is an exception for R^^^ ~ -^V^0^. Taking 
into account these canonical normahsations, the expansion of the action now takes the 
form: 

+ ■■■ + M^^rn-lH.^-2 i9'hr- + M2JH^rn-^ (9'hr-V^4>' + •.•), (23) 

and we see that the strong couphng scale, given by the lowest scale suppressing non- 
renormalisable operators, is = {MpH^y/^ . 

It is now straightforward to translate this strong coupling energy scale to the Vain- 
shtein radius, defined as the distance where the perturbative expansion of the spher- 
ically symmetric solution, corresponding to a central mass source, breaks down. In 
terms of Feynman diagrams the perturbative calculation of this (classical) solution 
amounts to the calculation of the tree level expectation values of r and 0, induced by 
the mass source with coupling 

where the last equation follows from (jlUUp in the limit of short (physical) distances f = 
ra{t) <^1/H from the source. The Feynman diagrams for the first order expectation 
values of and are given in fig. la. The definition of the canonical excitations is 
such that the only dimensionful quantity entering in the short distance propagators 
is the energy /momentum g; this gives a canonical dimension 1 for r'^ and a canonical 
dimension for (f)'^. With this in mind one can easily estimate the short distance 
behaviour coming from a certain diagram, by taking into account that the source has 
a dimensionless factor M/Mp for the r^-coupling and a dimension one factor HM/Mp 
for the ^'^-coupling. From the first order diagrams we then find 

rm Ml HM , , 

The dominant higher order corrections will come from Feynman diagrams con- 
taining only the third and fourth order vertex of (j23|l . since these are the only ones 
suppressed by A^. We depict some of these higher order graphs in fig. lb and fig.lc. 
Taking into account the relevant powers of A^ that come in front of the vertices, we 
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(a) 





(c) 



Figure 1: Tree level tadpole Feynman diagrams for the expectation values of t'^ and (j)'^, 
induced by a mass source: the first order diagrams (a), some leading higher order diagrams 
for the expectation value of t'^ (b) and a leading higher order diagram for the expectation 
value of (f)^ (c). The full lines represent propagators of the fields Xiji'^i and r. The 4) 
propagator is represented as a wavy line. 



now find that the higher order (m > 1) corrections go hke: 



r 



c(2m+l) 



My™ 1 1 



M 1 [GnM 



Mpf 
M2 



r4 



2m 



GnM 



2m-2 



(26) 



and we see that the expansions for and 



mV" 1 



Mr, 



diverge for distances f smaller than 



A. 



(27) 



Notice that we have not written down an expression for the even orders in the ex- 
pansion and of the odd orders in the 0"^ expansion. The reason is that the Feynman 
diagrams corresponding to these orders unavoidably contain one or more of the "sub- 
leading" higher order vertices in ()23|). suppressed by a mass scale bigger than A^. The 
coefficients of the even higher order terms for the expansion, for instance, will then 
go as 

.c(2m) ^, 21^2. / ^i^mi^ 1 ^28) 



r 



M 1 fGNMH\ 
{HrY ) 



Mpf 



where s < 3; and still show a convergent behaviour at distances r < ry, where the 
total series diverges. 
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4.2 Explicit check: the Vainshtein radius 

In the previous subsection we derived diagrammatically, for general Lagrangians of the 
form F[R, Q — 4P], the short distance behaviour of the perturbative expansion of the 
spherically symmetric solution corresponding to a central mass source. In terms of 
the original fields r and 0, our results (|^. (|^. for the expansions of and t'^ 
become: 

_ G^M {GMM)^^^{GMMf 

Here Xj are dimensionless coefficients and we have used = + (2 + 0{H'^f'^))T. 

Let us now consider a specific model, the action ((H) with n = 1 (and 6 = —4c), 
and explicitly check that the first orders in the expansions above do indeed exhibit this 
behaviour. At linear order the equations of motion, obtained from (f7H|) . (f7^ . (jHUI) and 
(I24j) have the form: 



16.G„ . V - 2.)J + ^( jg;. + . . .) j + . . . ^ (30) 

(f ) - (2§^<* - ^-)) - ..).... ^ -^^"W- (31) 

We will solve these equations in a perturbative expansion in rH . At short distances, 
the dominant terms in the equations are the V^-terms. These terms originate from the 
term oc i?*^^)^ in the quadratic Lagrangian. At first sight, this term would seem to give r 
a propagator oc However, as we discussed in the previous section, the propagator 
has in fact the canonical oc behaviour. In the equations of motion ()30I3H) this 
is refiected in the fact that the V^-terms are proportional to each other, which has 
consequences for the leading order (in fH) behaviour of the solutions. Looking at the 
2nd equation one would naively infer 

r, oc a{tf'^GNMH^r{l + 0{r^H^)) , (32) 

since V^r = ~'8)Tt5^{x). However, this is in confiict with the first equation. The proper 
ansatze will be: 

^(1) ^ a{tfn^^{Y,^^^^^^{f^H^r) (33) 

^ m=0 

0^^) = a{tr'^{Y.{dm + c^W'H'r). (34) 

m=0 
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The coefficients are written in this pecuhar form, because we then have: — 2r = 
^ c„(. . .). Plugging in this ansatze and solving the equations order by order we find: 

(...)V^J = ^ co = 0, (35) 
{...)S%x) = ^ do = ^ ,ci = -A_^, (36) 

(...)- = - di = -- = + (37) 

The dimensionless constants Ci and C2 are defined as 

Ci = 5 + ^ and C2 = ^ (38) 

mo mo 

where mo and 5 have been defined in the third section. So we find for the spherically 

symmetric solution, at first order in G^M: 

^(1) ^ ,(,)3/24^(2-(^ + ^)r-i7^ + 0(r-i7^)). (40) 

Notice that the leading behaviour is precisely the one we expect in scalar-tensor theories 
of gravity, but with a rescaled Planck mass — > CiM^. We have also calculated the 
subleading correction because this is required to obtain the leading correction for the 
second order 0^^^ and r^^^. These second order corrections are the solutions to the 
equations: 

(1) /XC\(2) 



(^|)'"(x«,0«) = -(g"(x«.,«). (42) 



An explicit calculation gives: 

(2) 

/ An T\/r\ ^ 



(1)'^' (.a.,,<0) ^ . (!^)\,)3.1iz^, , , . 



with: 



_ 732C2 + 79C| + 96Ci - 724C1C2 - 96Cf 
^ ~ 60C2(12 - 12Ci + C2) ^ ^ 

24OC2 + 29C| - 240Ci - 26OC2C1 + 240Cf 
' ~ 30C2(12 - 12Ci + C2) ■ ^ ^ 
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Notice that the leading terms in these equations have the same ratio as the V^-terms 
in the first order equations of motion. This can be traced back to the fact that they 
come from the third order vertex in (j^Hj) . And now, in contrast to the first order case, 
the leading behaviour will be solution of V^r"" = ^ a = 4. Writing then the 
expansion as 



(^) ' fiTP {t^^if'H'r^ , (47) 
_ „(,p('!G^y_J_/'g(;^ + ,^,,f.^.)™\ (48) 

^ ^ ^ \m=0 J 

and solving the equations order by order, we now find: 

(...)^ = ^ eo = (12Ci - 12 - C2)/(4C2) (49) 
(...)^ = ^ /o = eo, ei = .... (50) 
So the second order terms in the perturbative solutions are: 

rW « a(i)'/^^^^ (51) 
= 12Ci -12 - 2C. ^^^^3/. (^jGgjiy _L_(i ^ o(fV)) . (52) 

We will not obtain the actual 3rd order correction, but just show that its leading 
behaviour is: 0^^-* = 2t^^^ oc (GnM)^ / (r^H^). We find that the third order corrections 
are solutions of: 



a(t)3/2 (4^^^/)^ /12(8(Ci-l) + 3C2)(12 + C2-12Ci; 



IQttGn (3Ci)3fiiif6 ^ c'2(C2 + 4Ci - 4) 



> / 

a(t)3/2 (4G'^M)3 / 6(16(1 - d) - 3C2)(12 + C2 - 12Ci) ^' 
leTrGjv (3Ci)3fii/76 1^ ^'2(^2 + 4Ci - 4) 

Now the leading terms in the right hand side of the equations of motion do not have the 
special ratio that appeared in the calculation of the second order correction. This comes 
without surprise, since they originate from the fourth order vertex in that does not 
contain an R^^^ term. This means that analogously to the first order case, the leading 
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behaviour of t^^\ (p^^^ will be solution of: V^r" = r~^^ a = —9 and we indeed find the 
leading behaviour 0*^^-' = 2t^^^ oc (GnM)'^ / (r^H^). Notice that if we would have looked 
only at the second order correction we would have inferred a wrong Vainshtein radius 
(~ {GnM / H'^Y^'^). This was anticipated by our feynmandiagrammatic derivation: the 
even orders in the perturbation series for r (and 0) come with bigger suppression scales 
than the strong coupling scale A^, the scale that determines the true Vainshtein radius. 

5 Conclusions 

Models that involve inverse powers of the Kretschmann scalar Q in the action not only 

modify cosmology at late times, but also Newtonian gravity at long distances. The 

reason behind this behaviour is clear: in cosmology, when the Hubble scale is H ^ fi 

the extra term that appears in the Friedmann equation is suppressed by powers of ^/H 

and therefore negligible. Analogously, for spherically symmetric solutions, the extra 

term in the equations is now suppressed by powers of r/rc, and this suppresses the 

modification at short distances where the usual Einstein term dominates. Notice that 

the structure of the spherically symmetric solutions for theories of the type in a 

spacetime that has constant curvature asymptotically is such that the scalar curvature 

goes to zero at short distances and to a constant far from the source: this also reflects 

the complete shielding of the modification of gravity at short distances from sources 

that takes place in these theories. 

We have also analysed the particle spectrum of the theory. In this kind of theory, 

that involves up to fourth order derivatives of the metric on its equations of motion, 

one can expect generically eight propagating degrees of freedom in the gravitational 

sector: two for the massless graviton, one for a scalar excitation and five for a massive 

ghost. We have seen that the massive ghost is absent as long as the action depends 

on Q and P only through the combination Q — 4P ^. Moreover we have seen that the 

mass of the scalar field depends on the background curvature. This enables us to offer 

a physical picture in terms of particle excitations of the behaviour of the solutions: the 

extra scalar field becomes massive, and decouples, at short distances from sources. At 

long distances its mass decreases to a very small value ms ~ H and its effects can be 

^This is precisely the combination appearing for instance in the Gauss-Bonnet term. It is amusing 
to note that although the addition of this term to the action does not affect the equations of motion 
in 4D, the addition of its inverse yields a ghost free long distance modification of gravity. 
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noticed. So we are led to a situation in which gravity is well described by conventional 
General Relativity at short distances from sources, r -C Vc, or in regions where the 
curvature is much bigger than fi^. On the other hand, at distances r ^ ry, or in 
regions where the curvature is low, gravity is well described by a scalar-tensor theory, 
obtained from the linearisation of the action in vacuum. In both cases the force of 
gravity will exhibit the characteristic falloff, albeit with different proportionality 
constants. But there is still the intermediate range, Vc < r < ry, where there is a 
transition between the two regimes, and where we can expect deviations from the 
dependence in the force. To understand qualitatively this transition we can consider the 
modification to the Newtonian potential produced in scalar-tensor theories of gravity 
and argue as follows. A massive scalar with mass ms and gravitational couplings will 
typically contribute a term in the potential that goes like 

Vs ~ ^e--^ (55) 
r 

with the characteristic exponential Yukawa suppression of forces mediated by massive 
fields. If the mass is of order H we would expect that the exponential will be well 
approximated by one everywhere. But we have seen that in our case the situation is 
more complicated because the mass of the scalar (and the Planck mass) depends on 
the background curvature. As we have shown, in de Sitter space the mass of the scalar 
is given by 

p2n+2 

"I ~ ^ (56) 
where Rq is the background curvature (or, more precisely, the vacuum expectation value 
of a function of R, P and Q). We expect then that, locally, the mass of the scalar 
will have this dependence on the background curvature. If we now take into account 
again that in spherically symmetric solutions, a generic scalar involving n powers of 
the Kretschmann scalar is of order ~ (9 ((GM)^" /r^") at short distances, this 
translates into an r-dependent contribution to the mass of the scalar for actions of the 
type (P) that goes like 



^Source 



ms{r) ~ O ((G'M/r^)"^' . (57) 



Now we see that the scalar mass has a contribution that grows at short radius. We 
can obtain the distance as the distance at which the "exponential factor" in the 
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contribution of the scalar to the potential is of order one, i.e. 

exp[—ms{r)r] ~ 1 for r ~ r^. (58) 

For smaller distances the would-be Yukawa exponential factor suppresses the scalar 
contribution to the potential, effectively decoupling this degree of freedom. 

We have seen how this model shares many characteristic features with other theo- 
ries that also offer a long distance modification of gravity, like massive gravity or the 
DGP model. In all of them the perturbation theory valid at long distances from sources 
breaks down at short distances. This characteristic shielding behaviour makes it pos- 
sible to avoid the vDVZ discontinuity, and we recover an acceptable phenomenology 
even if there is an almost massless scalar field in the spectrum. In the models that we 
have studied here, however, we have a perfectly well defined theory above the cut-off 
of the linearisation, in contrast with massive gravity or other proposed long distance 
modifications of gravity. One might still worry about the possible appearance of more 
unwanted ghost-like degrees of freedom if one takes into account the higher order terms 
in the expansion in powers of the fiuctuations of the gravitational equations, as hap- 
pens in massive gravity |22]- However, one can expect that the mass of any such state, 
coming from the scale of suppression of non-renormalisable operators, will be at least of 
the order of the cut-off of the linearised theory, and therefore from an effective theory 
point of view this would not really be a problem''. As a final cautionary remark we 
will also mention that although we have proved the existence of acceptable solutions 
in the short distance, r <^ Tc, and long distance r ^ ry regimes, we have not proved 
that one can get a consistent matching between them through the "nonperturbative 
region" (similar potential problems have been pointed out for the case of massive grav- 
ity in [21). We regard this issue as more model dependent, and we hope to address 
the behaviour of the solutions in this intermediate regime more quantitatively in the 
future. 

^Similar remarks can be made regarding Einstein gravity. Quantum corrections generate higher 
order curvature terms in the action that inevitably produce ghosts upon linearisation. But the mass 
of these apparent states is typically of the order of the Planck mass, so they do not really represent a 
problem since one can not trust the linearised action in this regime. 
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Appendix A. The gauge invariant degrees of freedom 

In this Appendix we will offer the decomposition of the fluctuations of the metric 
in terms of gauge invariant combinations. For the actions we will be interested in, 
Minkowski spacetime will not be a solution so we have to expand on another metric. 
We will take a cosmological FRW background metric: 

ds^ = gl^x'^dx" = -dt^ + a{tf{dx^ + dy'^ + dz^) . (59) 

Writing g^^ = g^^^ + h^^ and leaving the background metric {g^^^,) fixed, an infinitesimal 
coordinate transformation x ^ x + gives: 

^U(^) = V(^) + v^e. + v.^M = V(^) + vje. + "^X + o{H) , (60) 

where the superscript indicates the covariant derivative on the background metric. 
We will use the same decomposition as in [22] to obtain the gauge invariant degrees of 
freedom contained in the perturbations: 

hm = V' 

hoi = Ui + diV (61) 
hij = Kij + {diSj + djSi) + didja + 6ijp . (62) 

Here transverse-traceless tensor; Vi and Si are transverse vectors, while the 

other fields are scalars under global 3D rotations. At lowest order the tensor modes 
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are gauge invariant, while the other modes transform in the following way: 

^ ^ V + 290^0 

Ui ^ Ui + ido-2H)Cf 

V ^ V + {do - 2H)7] + ^0 

Si ^ Si + ^J' (63) 

(7 — > 0" + 

where is the transverse part of ^j, rj the longitudinal part {di^i — V^rf) and H is the 
Hubble constant of the background metric. There is one gauge-invariant combination 
in the vector sector: 

w^ = u,^ (So - 2H)s, + 0{e) , (64) 

two in the scalar sector: 

= iP-2dov + do{do-2H)a + Oih^), 

r = p + a^H{2v-{do-2H)a) + 0{h'^), (65) 

and finally there is one gauge-invariant (transverse-traceless) tensor: 

Xij = Kij + 0{h^). (66) 

Obviously, these gauge invariant modes are the only ones that will show up in the La- 
grangian. Furthermore, the background is invariant under 3D rotations, which implies 
that at the quadratic level a generic action will have the following form: 

S^^^ = jd^x (x^jO^X^J + Wid^w, + + tO^t + 0630) , (67) 

where the operators O depend explicitly on time and contain time derivatives and 
three-dimensional Laplacians. 

Appendix B. The quadratic action 
Bl. The case of Einstein gravity 

As a warming up exercise we will briefly review here the situation in Einstein gravity 
for the gauge invariant decomposition we are considering. In this case F[R, P,Q\ — R 
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and the linearised action in flat space one gets after integrating by parts and discarding 
total derivatives is 

^^'^ = {-X^J^\i3 - '^w.V'w, + 40 VV + Grd^T - 2rVV) . (68) 

b47rGAr J 

where is the d'Alembertian and is the three-dimensional Laplacian. At first 
sight the action doesn't even look Lorentz invariant. But this is a constrained action 
and the propagating modes do have Lorentz invariant dispersion relations. The field 
does not represent any propagating degree of freedom since it appears in the action as 
a Lagrange multiplier. And it is clear from the equations of motion for Wi and that 
the vector and the scalar r also don't carry any propagating degrees of freedom: 

55(2) 



50 

55(2) 



oc V'r = (69) 



oc ^^Wi = 0. (70) 



6wi 

The only propagating modes are the two components of the massless graviton contained 
in Xij that have a Lorentz invariant dispersion relation. The solution to the equations 
of motion for and r coupled to sources do however generate Newton's potential. 

B2. The (dis)appearance of ghosts in the tensor/vector sector 

For Lagrangians of the form C = F[R, P, Q] we can expect up to fourth order derivatives 
in the action. So focusing for instance on the tensor modes we can expect, in the FRW 
background we are considering: 

jd^x XijOxXij = jd^x {x^,F^{t)Xij + X^o {Fl{t)V^ + F^{t)) x^, 

+X^, {F^m' + F^m' + F°(t)) x^^ , (Tl) 



where a dot denotes time differentiation. If is different from zero, the field Xij 
will contain 4 degrees of freedom, grouped into 2 degenerate pairs. From general 
considerations we know that the residues of the poles in the propagator will then come 
with opposite sign |7j, so one of the pairs will be a ghost. For -F2 ('^) find: 
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So we find that for actions of tlie form F[R, P, Q] = F[R, Q — AP] there is no ghost in 
general FRW backgrounds^. From now on we will consider that the action is indeed 
a function of P and Q only through the combination Q — AP. Also, we will restrict 
our analysis to de Sitter space in what follows since in a generic FRW background the 
results will be qualitatively similar, but the equations become more involved. So we 
take a{t) = c^* and we find useful to define the constants: 

Ci = 5 + ^ and C2 = ^ (73) 

ITLq ITLq 

where tjiq and 5 have been defined in the third section. After introducing some nor- 
malisation factors ^ : 

Xij a{tY^^Xij, 



and defining 



(75) 



a(t)2 ' 

we find that the operators in the tensor/vector sector of the quadratic action (j67|) read: 



= -fv^ (77) 

So we find that at the quadratic level, the excitations in the vector/tensor sector are 
the same as for ordinary gravity, but with a rescaled Planck mass: CiMp. As 

long as Ci > 0, the energy of the massless spin 2 graviton will be positive, or, in other 
words: the vector/tensor sector will be ghost free. We discuss the content of the scalar 
sector in the next section. 



^Notice that when this is the case we also ehminate two degrees of freedom in the vector sector 
and one in the scalar sector that together form the five degrees of freedom of the massive ghost. 

^When we write some of the excitations explicitly in section 4, we include these normalisation 
factors. 
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B3. The propagating degree of freedom in the scalar sector 

In the scalar sector, the operators appearing in the quadratic action are (again, assum- 
ing that F[R, P, Q] = F[R, Q - 4P], in a de Sitter background): 

-H (SC, + gc.) a„ + 3H^ (|C, + ^C,) , (78) 

^3 - -iS-t^-t^^-«'^(i^'+il^')- 

At first sight the scalar sector looks like a mess. Naively one would expect 3 propagating 
degrees of freedom. Also, as in the case of Einstein gravity, the action does not look 
Lorentz invariant. We will show explicitly now that there is in fact only one degree 
of freedom, and its dispersion relation is, as it should be, Lorentz invariant at high 
energies or short distances. 

To get the true degree of freedom, the easiest thing to do is to look at the equations 
of motion. We find that a suitable combination of the equations of motion for and r 
clearly exhibits a constraint of the theory: 

where O4 and O5 are defined as: 



Using this constraint in another combination of the equations of motion, namely 

5S^ ^ SS^^ 1 

6t ^ S(f) IUTTLjjv 



were now Oq and O7 are defined as: 



48 //M " V 4 G 



■'2 



- -i^ i He, + iv^ + 1 , (85) 
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one arrives at: 

^5o^-V^-i^^(| + 16|))r = 0. (86) 

So there is only one scalar degree of freedom, with mass m| = —H^ + 16^^ 

We will now determine the sign of the energy of the propagating mode. One way 
to do it would be by calculating the "energy-momentum tensor" of the gravitons, 
T^i, = — (l/87rG'Ar)G'|S, and plug in the solution. Here, G^S stands for the second 
order term in the expansion in powers of the fluctuations /i^j, of the full gravitational 
equations. We will follow a different strategy and look at the sign of the residue of 
the pole in the propagators. In particular we will be looking at the situation for large 
energy and 3-momenta: u!,k = k/a{t) ^ H, where the mode respects the Lorentz 
symmetry and the propagators are diagonal in Fourier space. We will also discuss the 
behaviour of the propagators for large off shell 4-momenta, since this information is 
useful for the discussion of the strong coupling scale in section 4. After a httle algebra 
we find the following expressions for the propagators: 



-1 



P<pT — ^ — -Ptt (87) 



Here the propagators are defined in such a way that for instance: 

< (t){x)T{y) >o= STrGjvi < x\P^r\y > , (89) 

where < (p{x)T{y) >q is defined as the tree level expectation value of 4){x)T{y) one 
obtains from the path integral. Again, at a first glance the propagators seem to contain 
more than one degree of freedom. Naively, we would think that for high energies, we 
can forget about the explicit time dependence and ignore the fact that the operators do 
not commute ([9o, V^] ^ 0). For the rr propagator we then expect to find in Fourier 
space: 



Prr{'^,k') 



02(0;, k^)Os{u;, k^) - 101(0;, k-^yO^iuj, k^) (. . .)^' + (■ ■ O^' + (■ • 0^' + (■■■)' 

(90) 
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implying 3 propagating modes. In reality, the coefficients in front of uj^ and uj'^ cancel 
out, but to obtain this second cancellation one has to take into account the time 
dependence of the background reflected in the non-commutativity of the operators. 
The easiest way to do this is by using some identities that follow from equating the 
operators in front of r and (j) at both sides of equations (jHT|) and (jHH)- Using 

_ 1 
5t \Q-kGn 
6S^^^ 1 
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N 



2O2r + Of0j , (91) 
(2030 + Oir), (92) 



we then obtain from 



202 + 6501 = 66 + 67, (93) 



0^ + 20^0^ = O 



7 ■ 



And from (jHTj) we get: 



1 + 0405)01 + 20402 = -SH^Ci, (94) 



2(1 + 0405)03 + 0401' = -m'Ci. 

Some manipulations of these equations learn us now that: 

1 At A_i a 3 



62 - -6^63-161 = --C,H^ (64 163 16466) . (95) 
So the rr propagator can be written as: 

where the last equation holds for high energies/momenta uj,k ^ H. Just as in the 
case of the spin 2 graviton, we see that the scalar will have positive energy (or positive 
norm), for Oi > 0. 

For feymandiagrammatical calculations, it is useful to project out the propagating 
mode and work with instead of 0, where we define the former as 

= 0+^r. (97) 

The 00 propagator then reads: 

P^4> = ■ (98) 
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and the (pr propagator is zero. Notice, that for high energies/momenta the 00 propa- 
gator has the non-canonical behavior 

^^0---^^, (99) 

which is subleading with respect to the canonical rr propagator. When dealing with 
the issue of the strong coupling scale in section 4 we need the explicit high energy 
behaviour of O^^Oi. 

O-^Oi ^-2(2- 3V-^(9^ - 9HV-^do + 9HV-^d^^ . (100) 

Notice that, on shell {uj^ = k"^), we have 03^^01 ^ 2, or = — r + 0, consistent with 
the equation of motion r -|- = for the propagating mode. 
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